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The complex time WKB approximation is an effective tool in studying particle 
production in curved spacetime. We use it in this work to understand the formation 
■ of classical condensate in expanding de Sitter spacetime. The CWKB leads to the 

emergence of thermal spectrum that depends crucially on horizons (as in de Sitter 
spacetime) or observer dependent horizons ( as in Rindler spacetime) . A connection is 
CO ' sought between the horizon and the formation of classical condensate. We concentrate 

^1 on de Sitter spacetime and study the cosmological perturbation of A; = mode with 

various values of m/Ho. We find that for a minimally coupled free scalar field for 
m /Hi < 2, the one-mode occupation number grows more than unity soon after the 
physical wavelength of the mode crosses the Hubble radius and soon after diverges 
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O ; as N{t) ~ 0(l)[Aj,h^,(t)/J/(7^]V"''-i/4^ ^jjgj.g ^ = (9/4 - mV-ffo)^''^ The results 

CNj , substantiates the previous works in this direction. We also find the correct oscillation 

' and behaviour of N{z) at small z from a single expression using CWKB approximation 

^ ' for various values of m/Ho- We also discuss decoherence in relation to the formation 

of classical condensate. We also find that the squeezed state formalism and CWKB 
^ ' method give identical results. 



X 



I. INTRODUCTION 



^ . An important aspect of cosmological perturbation is that they are of quantum origin but eventually 

they become classical as they are (or source) responsible for classical density perturbation. A general 
trend to approach this problem is to consider the interaction between 'subsystem' and 'environment' 
such that the subsystem decoheres during evolution. The quantum fluctuations that originate in 
the early universe thereby undergo quantum to classical transition being stretched to astronomical 
size due to inflation. The works in this direction will be found in [1-14]. The decoherence due to 
particle production is another approach in which it is argued that the particle production during 
evolution decoheres the quantum fluctuations. Here we need not require to introduce interaction 
between the subsystem and the environment. Calzetta and Hu [15] studied decoherence of a mean 
field by its quantum fluctuations, and Calzetta and Gonorazky [16] considered fluctuations in a model 
with nonlinear coupling. In [17] decoherence is studied through particle production. The initial works 
in this direction are due to Guth and Pi [18], Lyth [19], Polarski and Starobinsky [20], and Lesgourges 
et al [21]. 

Let us try to understand how decoherence is obtained through particle production. Guth and Pi 
showed that at late times the wavefunction of long-wavelength modes behave as classical Gaussian 
probability distribution whereas Lyth pointed out that the wavepacket of massless modes in de Sitter 
space does not spread and this property is attributed to classicalization of fluctuation. In [20,21] one 
looks at the free fields and the classicalization of the modes occur due to particular property of de Sitter 
spacetime. In this approach the decoherence emerges without invoking any interaction which is deemed 



1 



essential in usual decoherence studies. One easiest way to study classical nature of a quantum variable 
is to look at the condition under which it becomes classical. For example, the canonical coordinate 
will behave like a classical variable if it approximately commutes with the canonically conjugate 
momentum. In expanding spacetime we are to choose a variable whose emergence is due solely to 
quantum behaviour of the fluctuation and in someway is related to classical density perturbation. One 
such quantity is particle occupation number. Assuming that vacuum fluctuations (that are supposed 
to act as seeds of density pcrturbation)in de Sitter spacetime during evolution produce particles, as is 
usually the case in expanding spacetime, we can search for decoherence looking at how this occupation 
number behaves at late times. If the number of produced particle is enormously large, we can argue 
that the situation is ripe enough for the formation of classical condensate. In other words decoherence 
is accomplished through particle production. This approach is somc^tinics known as 'decoherence 
without decoherence'. In de Sitter spacetime the classical nature of fluctuation is usually 'assumed, or 
argued for in passing, by invoking Gibbons-Hawking radiation, or the horizons in de Sitter space'. It 
is worthy to mention that a system with zero dissipation may well be in a pure quantum mechanical 
state, even though the occupation number is large. For this reason one may have some apathy to 
the idea of decoherence without decoherence. However the present attempt may be continued in the 
manner of [22] to the standard picture of decoherence in which the large fluctuations plus the inflaton 
zero mode assemble into a new effective field. This effective field behaves classically and it is this 
object that rolls down. In this sense the study of formation of classical condensate might be a viable 
attempt to study decoherence. 

We, in this work, will follow [17] but will adopt the CWKB (complex trajectory WKB) method of 
particle production. The reasons for doing so will be clarified in the text. We consider a massive field in 
de Sitter background and study the conditions under which the vacuum fiuctuations in this spacetime 
evolve into classical fluctuations. We also discuss how the results obtained in the present work can be 
related to the standard method of studying decoherence through squeezed state formalism. 

To start with we discuss in section II the field theoretic approach to calculate the number of 
particles produced in a A; = expanding de Sitter spacetime with a Hubble parameter Hq for various 
values of mass of a free, scalar field. In section HI we obtain the CWKB extension of Nk that we adopt 
for numerical calculation. In section IV we discuss the numerical results obtained through CWKB 
analysis of section III. In section V we end up with a conclusion where we discuss the results obtained 
in the present working relation to squeeze state approach [4,24]. 



II. DE SITTER EXPANSION 

Let (f) denotes a scalar field of mass m, and S{r]) the scale factor for the Robertson- Walker spacetime 
with conformal line element. For the de Sitter case S{r]) = —j^- The Lagrangian in a conformal 
spacetime reads (minimal coupling) 

^/^L = ^5V(5,.</')(a.<^)-^54mV (1) 
with the metric, — S'^rj^i,, rj^i, = diag(l, —1, —1,-1). The field equation now reads 

^r^^''d^{S^d.cj^)+S^m'cj> = Q (2) 
We now decompose the real scalar field into the modes as 



Sj (27r)3 

The mode function now satisfies the equation 



e-'=-K^(,?)ag + <g(,?)4] 
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and the Wronskian 



For de Sitter spacetime, we get 



The Hamiltonian corresponding to (1) is 

which in terms of creation and annihilation operators and the mode functions reads 
where in de Sitter time we find using cqn (4) 



^niv) =\^-+ drij uj:{v) + a, j «i(77) + cji{v)u-,{v)ul{v) 



with 
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The conformal time Heisenberg operators ag, at in (8) can be written equivalently in terms of the 
77 dependent creation and annihilation operators 6^(77), ^^^(77) that makes the Hamiltonian diagonal. 
Introducing the mode function Vjg(?7) with the field decomposition 
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we get for the Hamiltonian 



H ^ - 
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(12) 



(13) 



For E'-^irf), F^{v) '^^ g^t similar expressions as in (9) with Uj^ replaced by Vj^. The time dependent 
Bogolubov transformations that diagonalize the Hamiltonian i.e., Fj^iv) = ^ow satisfy the relations 



<gW = PUV)V^{V) + aUv)v*_^{v) (14) 

and 

bl^iv) = + (^Uv)al^ (15) 

with |afe(ry)p — \Pk{v)\'^ = 1- The reverse transformation of (14) and (15) can now easily be found out. 
One can easily check that the new creation and annihilation operators satisfy the standard algebra, 
and they act on a time dependent vacuum as 

fcfe(^)|0>r,-0, bl{f^)\0 \k >„ (16) 

For diagonalization, we see from (13) that, the mode function Vj^{r]) has to satisfy the condition 
Using (9) with Uk replaced by Vk we get from the above relation (17) 

Using this result in the expression of E'Jrf) we get 



(18) 
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Using (8) and (13) we now get 

<0|26t(,)6,^(,) + l|0>=|||| 

Hence 

N,.(n\ = 

2 EL{^) 2 

where 



k^ 

E'Jr,)=2<^l{n)\vA-' (19) 



(20) 



mv) (21) 



Nk{ri) =< hib-k >= mv)f (22) 

is the number of particles referred to the initial Fock vacuum of the modes ug. The result (20) can 
also be deduced from F;.{r]) = F^{ri){al{r,)f - 2Ek{v)al{v)(3k{v) + F^{v){(3k{v)? = 0- This gives 

ak{ri) = PI ( ^1+"'= ) . Using la^p - |/3fep = 1, we get back the result (20). Further using ^^(77) = 
and E'f,{ri) = <jJk{'n)^ we will find \v^{ri)\ = (2ujki'r]))~ and hence (19) gives E'f,{ri) = oJkif])- It should 
be pointed out that the conditions ^^^(j?) = 0, E'j^{r]) = ujk{r}) do not determine the phase of Vk{rj). 
Using (9) and (18) we now get 

Nk{Tl) = -^^^ ^''^"fe^i + ^"^^i + ^kiv)UkUl _ 1 



4 ^liv)\v-k\^ 2 

1 + \"'f + ^i'^>i"i;""i: + '^n"",:"k) + 7Fl"fcl" 1 

4 ^l{v)\vj;\' 2 



(23) 
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'^fc = '^fc - Z2 = + ^ (25) 



So far the calculation is exact. Now to calculate Nk{rj) wc are to find ut from (6) and Vk from (18). 
We note that in (23) Wk is always real and positive but in the mode eqn (6), we have the oscillator 
like equation 

Uk " + ^luk = (24) 

where 

2 _ 2 a - 2 
— 5- — K H 5 

with a = /Hq. Thus for a> 2, uJk is everywhere real but for a <2, ujf. turns out imaginary when 
77 — > 0. In the literature most of the workers deal with real frequency but we will shortly find that 
the emergence of complex frequency has interesting features when one calculates Nu through CWKB 
technique keeping the content of (23) unaltered. In our earlier works [25 29] wc studied particle 
production through CWKB trajectory analysis with successful results. We adopt the same approach 
here. In Mijic's work, in defining the annihilation and creation operators as in standard harmonic 
oscillator prescription, he used \LOk\ instead of Wfc. This allows in dNk/drj to have a part proportional 
to CjI — pfep. Some argue this to be an improper way to define the operators. In this work we 
eliminate this objection. Had one started with the standard harmonic oscillator like definition for 
a and a\ one would not get the part that we just mentioned. We have tried to eliminate all these 
objections and explicitly show the role of complex trajectory through the use of WKB mode solutions. 
Moreover, in [17] Mijic considered superposition of two WKB like solutions to get U2 in evaluating 
Nk at z ^ in contrary to the essence contained in (23). We give a clear physical arguments of this 
type of superposition in terms CWKB trajectory. 

III. COMPLEX TRAJECTORY AND WKB APPROXIMATION 

The use of complex WKB trajectory to calculate particle production is well known [25-30]. The 
method has been used to calculate particle production in de sitter spacetime [26] obtaining correct de 
Sitter temperature kT = Ho/2n. Subsequently the method has been used to study boson and fermion 
production using oscillating inflaton [31] with remarkable success that also settled some variant results 
in the literature. The reader may go through [27] to understand the basics of CWKB. In the present 
work we will quote some results of CWKB relevant to the context of the present paper. Before we pass 
to the CWKB calculation let us discuss Mijic's result in this regard. Considering (24) as a oscillator 
equation with time dependent frequency, and introducing 



Ofc = (-)n/ + ^^\Pk (26) 



4 = (H-)V^"^ + ^^^ (27) 

and with the number operator =< aj. > it is shown that 

^ j^i^(42 _ _ 1 M(42 + al) (28) 

In obtaining (28) one has to use Pk' = u'l = —w'f. Uk and replace Uk and pk using (26) and (27) in 



terms of ak and a\. 



Now it is evident from the expression (28) that we have two sources of particle production. The 
first term contributes when the frequency is imaginary and vanishes for real frequency. The second 
term contributes for all frequencies. Most of the authors consider the second term in evaluating the 
particle production while considering the evolution of in vacuum to out vacuum through Bogolubov 
transformation technique between in and out vacuum. With (omitting the subscript k on u) 

ui,2(7?) = , ! ^ e^'J '^'•f")'^'' (29) 
the number of particles produced is obtained as 

mv) = = ^-Vt'^'^i' + - 1 (30) 

2\ujk[i])\ 2 2 

Referred to (26) and (27), as noted earlier, some have apathy in writing down the creation and 
annihilation operators in terms of \LOk\ for obvious reasons. Secondly, the solution U2 is obtained in 
terms of Hankel function so that in the small r] approximation, one has to put some restriction on a 
values. Mijic obtained large particle production when a < 2 and oscillating for a > 2. Leaving 
aside some minor technical points, wc clarify the first point through field theoretic calculation as in 

(23) . We then use WKB modes to calculate U2 without using the exact solution of temporal equation 

(24) . This is the advantage of CWKB over any other method. Moreover the definition of vacuum in 
curved spacetime is very delicate. The WKB like mode solution serves as a good template for defining 
the vacuum in curved spacetime. As we require only U2 in calculating Nk, henceforth we call this 
solution as Uk in the following. 

We now take the WKB modes as 

u,{v) = -^^=e-'I-~''(^^^^ (31) 

as solution of (24). Now 

Uk^iUkil-2/r^^Luiy/^ (32) 

As iVkiv) is always real, we thus get 

ojUv) = \wk{v)f ^ \u>k\\l + 2/rful) (33) 

Using this in (23) we get k'^rf + a> 2. However this is not so important, it is mentioned here to show 
that is positive. Since 

\ ■ - Co' Cj' 

Uk = ^= ^ '^''^'^ , dr/Uk = -T^Uk - iCokUk i-e- {dr,UkUl + Ukdr^ul) = --^\ukf 
V2w ^Wfc uj 

From the definition of u)k we get wjj. = — thus, 

{driUkUl + UkdrjUl) = -^^\Uk\^ 

We can write 

|2 |2 L/ |2 L,/ |2 



„2 



'^k \'^k\ '^fe 

U,' |2 



2 



4^2 



\4\ 



For small 7] and neglecting the last term we have 



Thus 

I I I ] I / '^1 1 I 1 2 I 1 I 12 
*(,) = j " |,,p • 2 l^") 

At small T], the last three terms in the numerator cancel out almost exactly. For large t] the contribution 
of last three terms also can be neglected when compared with the first two terms. We fix the initial 
vacuum at large rji by taking 

^;fe(r?) = -=L— e^^-^f"^)" (35) 

Thus we get 
where 

IF, (r,\\ - - |(l+(a-2)/(fc,).)^)^/^| 

\^WU\ - - |(l+(a-2)/(fer,)2)i/2| 

We note some characteristic features of our expression (36). 

(i) It is well known that WKB approximation is valid for large 77. We show that WKB approximation 
as carried out in this paper, works very well in small 77 region. In other words our expression (36) is 
valid for small and also for large t]. So it is worthwhile to study these regions analytically. For large 
r], the prefactor |-Ffe(r?)| 1 since the vacuum is defined at rji — > —00. We then get back the result 
(30). 

(ii) In WKB modes to define the vacuua are well known, usually known as Parker's definition. These 
modes serve as good template for particle and antiparticle states. Not only that, our definition used 
in the present work is also being currently used in the literature [23]. However this type of WKB 
definition has also relevance when one discusses squeezing to understand decoherence mechanism. 

(iii) At 7] ^ 0, the WKB mode solution is not valid. Here there is a standard analytic continuation 
to keep the WKB mode even at 0. We first consider m = case in which analytical results are 
available for comparison. We get going to the variable z = —kr], the equation (6)as 

u'l{z)+u>l{z)uk = (37) 



and Nk{z) = ^\Fk{z)\ 



\^k{z)\\uk\^ + r.-j-r-Mz)\' 



(38) 



where now the prime denotes the differentiation with respect to z. Introducing a parameter A such 
that X^P = /3 we write (37) as 

u'az)-X'^u, = (39) 
where A is a large parameter and P{z) = 2 — z^ . The WKB solutions at ^ — > is now given by 
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Uk ~ e^"^ ~ z^^^ (40) 

where 

S=l{2-zy— (41) 



The validity of WKB solution (40) can now be verified through the exact solutions of (39). Doing 
exact solution to (39) we will find 

This gives us (40), keeping WKB mode solution even at ^ 0. With this change and evaluating 

S{z) = LJk dz 



we find after setting X P = P 



Nk{z ^ 0) ~ 2^ - 4 (43) 



where the extra factor z comes from the Fk{z) term. Without the \Fi.{z)\ term, wo will get Nk{z) ^ 
1/^2V2. If we do the exact calculations finding Uk{z) in terms of Hankel functions we will get eqn 
(43). This is a very important result. 

So far we have nothing remarkable new. If particle production in de Sitter spacetime has anything 
to do with decoherence, we must be able to reproduce thermal spectrum at large z and and large Nk at 
small z. We know that the existence of horizon associates with the de Sitter spacetime a temperature, 
which fact in our approach, will be taken into account through the introduction of complex trajectory. 
For w? /2Hq « 1, well satisfied for vacuum fluctuation, we have negative u"^ at small r]. Thus we 
have both growing and decaying modes, and if growing mode survives we have large Nk and this 
allows us to connect A^^ to the existence of large squeeze parameter, a condition for having effective 
decoherence through particle production [24]. In both the regions {m? > 2Hq and m? < 2Hq) we will 
find that the complex trajectory plays the crucial role. Let us first discuss the emergence of thermal 
spectrum. 

The emergence of thermal spectrum is straightforward in CWKB approach. We now write the 
expression of A'"^ in terms of z. We have 



Nk{z) = \Fkiz) 



\uk {z)\ + — ^ — \Uk{z)\ 



2|wfe(z)| 



We consider the case a » 2. In that case we have no real turning point. The turning points are now 
complex, given by 2:1,2 = i^s^- Keeping aside the pre-exponential factor in the WKB mode solution, 
we have 

Uk^e+'^''^'"'^'^ (45) 
u'f. = iujk{z)uk (46) 



Now time flows from large z to small z. Let us see how is modified when we take the contribution 
of complex turning point. Neglecting the vacuum contribution we have approximately (neglecting the 
variation of ijjh{z)) using eqn (44) 

Ar^(^)~ i|e'^A"'=(")''^|2 (47) 

Here Zi is an arbitrary initial real point from where the trajectory starts and ends at the real point z 
such that z > Zi i.e., we are moving backward in time. Now the trajectory starts from Zi, comes to 
a point zq (real) at which the particle becomes nonrelativistic, proceeds towards the complex turning 
point zi, being reflected again comes to on the real axis and finally goes to z. It should be noted 
here that we are considering Ufe ~ exp{iS) where S is the WKB action (not to be confused with the 
scale factor). We write this contribution in steps 

^ ^iS{zo,Zi)+2iS{zi,zo)+iS{z,zo) _ ^iS{z,Zi)+2iS{zi,zo) (4:8) 

with S{zf, Zi) being the WKB action between Zi to Zf. The integral has been evaluated in [26] to find 

,T / X -4Im f "^2(1+ » )i/2^ -4x^ /Atw 

so that 

Nk{z) ~ e-T (50) 

where T = Ho/2Tr with Boltzmann constant taken as unity. Thus we get the thcrmalized spectrum. 
The point zq on the real axis, at which the trajectory starts towards the reflection, can be found as 
follows. Starting with the scale factor S{t) — HQ^exp{Hot) and setting S{t) = ±«^, we find the 
turning points ti^2 = ^ijTi^ +^o'^'"'^^- Reti,2 we get a{Reti^2) = ^ = a{Rer]), where particle 
becomes nonrelativistic. We will discuss this case again when we consider the case » 2Hq 
separately. One important feature is that the expression of Nk{z) is nowhere is dependent on initial 
Zi where the vacuum is defined. This is sometimes referred to high squeezing in momentum i.e., all 
differences in the initial state are damped out by inflation. There is a general result to evaluate the 
thermal spectrum for any expanding FRW spacetime. The result is [29] 

/ Wfcda; ~ — / u)kda = iTrm/Ho (51) 
Jzo ida/dz),Jao 

The temperature comes half of that obtained through the exact treatment. Though all FRW expanding 
spacetimes show thermal spectrum, large particle production is not found in all such cases except for 
the de Sitter case. Note that in evaluating (48) we consider only reflected trajectory contribution but 
there may be another direct trajectory from Zi to z to interfere with reflected part and will generate 
oscillation in N^. We will discuss this separately. 

We will now discuss the CWKB results of particle production both analytically and through nu- 
merical calculation. We alrc;a(ly noted that the expression of Nk{z) does not depend on the initial 
moment Zi but it correctly reproduces the vacuum condition Nk{oo) = in the limit z ^ oo. 

The regime of interest is z ^ 0, when there is large expansion and the wavelength of the modes 
exceeds the Hubble radius. Here we are to consider two cases 

(i) >> 2Hq, real frequency with bounded classical motion and 

(ii) << 2Hq, imaginary frequency and rolling down the hill. 
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In the case of > 2Hq the frequency is everywhere real and there is no turning point for real z. 
However there are complex turning points at z = zLiM with M = y' w? / Hq — 2. It will generate a 
thermal spectrum with N ~ exp{—m/T) with T = Hq/2-k as already shown in eqn (50). However for 
quantum fluctuation, we are interested in those m— values for which frequency turns imaginary for 
/Hq « 2. In this case the frequency turns imaginary at zq = «"\/2\/l — "ni? /{2H'^). If 770 denotes 
the time when frequency becomes zero and r}c the time at which the fluctuation crosses the Hubble 
radius, we get 

S{m) = S{ri,)(2-'^\ (52) 



This relation implies that for the scalar field (supposed to have low mass << H"^) that drives 
the inflation, the one mode frequency turns imaginary before the mode crosses the Hubble radius. 
For Hq < rm? < 2Hq, the "rolling down the hill" occurs even at later times, well after the physical 
wavelength of the mode crosses the Hubble radius. In other words, it is the evolution for 2: < 1 that 
decides the fate of the fluctuation. We will discuss the case > 2Hq later as it involves the inclusion 
of complex trajectory. 



Case : < 2H^ 

For the case rn? < 2Hq we are to evaluate N{z) within the real turning points ztzg. So it is the 
evolution for z << zq that decides the fate of the fluctuation. Here again we will apply the CWKB 
method to keep the particle-antiparticle definition fixed with respect to initial vacuum. From (38) we 
get the expression 



Nk{z) = ^{e^'"'^ '"'^'^'''\Fk{z)\-l) 

We now evaluate the integral 

/(z) = y" Uk{z)dz = Mlog + M^) - - Mlog z + \J z^ + M2 

where = (2 — ^). With w? /Hq « 2, taking the imaginary part of I we finally get 



(53) 



(54) 



\Fk{z)\ 



2a + 4 



(55) 



where a = rrP' /H^. It is evident from the expression that at z ^ 00, N{oo) = satisfying the initial 
vacuum condition. In evaluating the above expression of Nq[z) we have not taken the variation of 
frequency i.e., we have taken 



(56) 



This No{z) is plotted in fig. 3 for various values of m/Ho and the production of large number of 
particles is quite evident. Let us see how this iV(z) is modified when we take particle production due 
to change in frequency. We now get 



Nk{z) = N{z) = ^\Fkiz)\ 



z'^ -2a + 4 



7(2^' 



1 + 



'|wfc(2)|^ 



(57) 
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This Nk{z) is plotted in fig. 4. This fig. 4 is in in conformity with the variation of frequency shown 
in the fig. 2. We obtain the fig. 2 using the expression 

kLi _ a-2 

where a = 
Case: > 2H^ 

If we consider a > 2, the frequency is real for all real values of z, so from eqn (36) 

\Ukf ~ (58) 

and hence in this case the particle production is entirely due to rapid change of frequency at small z. If 
we do not consider the turning points, which will be complex, the contribution like the first bracketed 
term in (57) will be now absent and hence we are to consider CWKB trajectory contribution more 
carefully. The importance of our approach is that the single expression for N{z) represents all the 
characteristic features of particle production using only WKB approximation contrary to various 
approximation as is carried out in [17]. 

In this case the integral of C0k{z) is always real and the particle number produced is much smaller 
compared to case < 2Hq. If we take the contribution of the complex turning point at ±iM with 
M = — 2 we find that the Nk{z) so reproduced does not result in any of the graphs shown in 

the figure (see [17], fig. 3). They obtained the graph with the Hankcl function approximated to small 
z values. Let us see what happens when we take complex trajectory approach. In this case there is 
no turning point at real z. The contribution of Uk may be considered as 

«.(.) = -^=l=e'^^-(^)'^^ (59) 

where Zi is the point where the initial vacuum is defined. The physical picture behind the CWKB 
particle production is Klein-paradox like situation , not in space but in time. Though this has been 
discussed elsewhere in details, we briefly discuss the physical picture to consider reflection in time as 
a process of particle production from vacuum. Consider the pair production from a time dependent 
potential V{t) as 

Vit) (e+ T)+ + (e- T)+ 

Here the uparrow denotes that both particles move forward in time with positive energy. Now ac- 
cording to Feynmann-Stuckleburg prescription a positive energy particle moving forward in time is 
equivalent to negative energy antiparticle moving backward in time. Thus the pair production process 
shown above can be written as 

V{t) (e+ T)+ + (e+ i)^ 

This equation implies that a reflection in time can be viewed as pair production [25]. This corresponds 
to a situation where referred to eqn (24), particle production occurs at 0Jk{z) = The turning points 
may be real or complex. For real turning points we have situation just like tunneling through a 
potential barrier (case m? < 2Hq); for complex turning points, the situation we encounter is like 
over the barrier reflection in potential problems (case > 2Hq). The reader may consult [27,31] 



to view the CWKB trajectories in tunneling and over the barrier situation in potential problems. In 
CWKB we consider (24) like Schrodinger equation not in space but in time. The wave function at 
any real point is written as contribution of WKB trajectories. In CWKB, the variable (here z) is 
treated as complex. The wavefunction at any real point can now be written as sum of real trajectory 
and complex trajectory contributions. The incident, reflected and transmitted trajectories are then 
considered in the standard way. For details the readers may consTilt [27]. In dc Sitter case, we have 
z as the time variable. Thus z — » — oo correspond to the situation where there is no particle and 
z — > +00 correspond to the situation where we have pair production. The turning points in between 
the two regions given by w{z) = specify the barrier region. Here exp[i S{z, Zi)] acts as incident wave 
and as it emerges from the complex turning point region, the wave is written as in eqn (66). 

Previously while discussing the thermal spectrum we consider only the contribution of branch point 
of 0Jk{z) and noted that the result is independent of Re t on the real axis from where the trajectory 
starts to get reflected. It is therefore necessary to consider the trajectories that go direct from Zi to 
z and also the trajectories that suffer reflection at the turning point and go to z. The trajectory that 
goes to Zi ^ z without any contribution from the turning point is 

The other trajectory starting from zt goes to the turning point zi and turning back from zi again 
comes to z. This trajectory contribution is 

Here {—i) has been introduced for reflection at zi and S = Cj{z)dz = real. We now introduce 
repeated reflections between the two turning points. Eqn. (59) is modified as [27] 

ul = Uk{z){l-i\Re\e-^'^) (60) 

where 

+ exp{-2TTM)) 

where the numerator comes from the repeated reflections between the turning points. Hence we get 

\u^^\z)\^ = |«fc(^)|2 (1 + - 2\R,\ sin 2^) (62) 
Now we find oscillation in Nj~[z). For a > 2, wc have now the result (since 6 is real) 
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(l + |i?c|^-2|i?e|sin 2^) - i (63) 



In evaluating (63), we have not used the variation of 5 with respect to z. If this variation is included, 
(63) is modified as 

Nk{z) = N{z) = i(l + - 2i?cSin 2(5) + — (1 - i?csin 2(5) - ujR^cos 2 5)]\'^ 

4 4 2 

1 uj' 

+ |[-— i?cCOS 2(5 + w + wiicSin 2(5]p (64) 

2 u) 

We plot this A''(^) in fig. 5 for various values of a > 2. Interestingly, (64) is of the form 
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N{z) = Aq{z) + cos 2 ^ + As sin 2 5 



(65) 



where Ao = -i + i(l + R'){2 + 

A — I ^ R 
- + 2 

Eqn (65) is of the same form as in [17]. However our coefficients are different because in that work z 
dependence in 6 is introduced in a different way. However we need not require any approximation for 
small z. The importance of our approach is that we need not require to evaluate the mode solutions 

and hence can be generalized to more complex cases. The present work is a good example in the 
cosmology where the CWKB paths play a definite role with a transparent physics. It should be noted 
that we worked always within WKB approximation and do not require to find the exact solution of 
(24) in terms of Hankel functions. 
In CWKB [27] U2 for > is written as 

U2 = e*^(^'^*) - ii?e-'^(^'^') (66) 

where e*'^^^'^') is the incident trajectory and Q-^S(z,zi) ^-^c reflected trajectory. Here we have not 
included the WKB pre-exponential factor for convenience. For w?/Hq < 2, the reflected trajectory 
contribution is very negligible, e*'^^^'^*^ grows due to imaginary Cbkiz)- For iti/Hq > 2 it is easy to 
show from the expression (66) that (62) follows with R = Rc with the vacuum chosen as 



Thus we see that the CWKB approach explains very successfully the cases a > 2 and a < 2. 



IV. DISCUSSION OF NUMERICAL RESULTS 

The numerical results are plotted in figs.l - 5. In fig. 1 we plot frequency square for various values of 

the parameter a = w? / H^. It shows that for a> 2 the frequency is unity at early times and changes 
rapidly as soon as the physical wavelength of the mode crosses the Hubble radius. The oscillator 
remains upright all the times. For a < 2 the oscillator turns to an upside down in the same regime. 
We note from (13) that there will be a term to particle production proportional to We plot it 
in fig. 2 in order to compare N{z) for the cases oJfc' = and uJk 7^ 0. To compare the contribution of 
changing frequency we plot in fig. 3 N{z) with u'^, = itUki^k for various values of a = jp- < 2. With 
the same value of a we plot N{z) in fig. 4 taking into account the frequency change while evaluating 

For a <2, the frequency turns imaginary for z « 1. Here the CWKB paths take crucial role. 
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FIG. 1. The evolution of for various values of a = /Hq. For a greater or less than 2 diverges 
after the physical wavelength of the mode crosses the Hubble radius. In the regime of large acceleration 
< Iferjl < 1, the effect of changing frequency will be important. 
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FIG. 2. Variation of ^ with kr). The peaks in the figure come from the zero of ui^ . For all values of mass 
the value of 7- exceeds unity shortly after the physical wavelength crosses the Hubble radius. 
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FIG. 3. The evolution of one mode occupation number in de Sitter space for different values of mass when 
the effect of — is neglected. 
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FIG. 4. The evolution of one mode occupation number in de Sitter space with the effect of ^ included for 
m? /Hq < 2. The figure is obtained from eqn (57). The peak at 2: = \/2 for a = is due to passing through 
zero as the oscillator turns into a upside down. The corresponding divergence for a = — 4 is not shown in the 
figure as it is outside the range. For a = 1.25, we also get the peak at the correct position using the same eqn 
(57). The effect of ^ is now clear when we compare fig. 4 with fig. 2 and fig. 3. 



Case a < 2 : Here we take a = —4,1.25,1.1,0 and plot the results in fig. 3 and fig. 4. The 
behaviour of N{z) for the cases a = —4, 1.1, is the same as in [17], the peaks also comes in the right 
position in conformity with fig. 1. In our case the magnitude of N{z) is larger in all the three cases 
and is obtained with CWKB approximation unlike the results in [17]. This lends support to the fact 
our result for the case a = 1.25 i.e., v = 1 is quite correct. Mijic [17] treats this case as special case. 
The peak at z = 0.85 corresponds to the zero of cj^ at this point. The difference that we observe may 
be ascribed to the following facts. When we calculate \(3k\^ using CWKB we get the last term of N{z) 
as in (30). Obviously the approximation made in [17] for Uk will not lead to a factor —i resulting in 
the factor ^. With our WKB choice of modes we get the factor —i i.e., we measure a with respect to 
the in vacimm with zero particle content i.e.. a;| >= 0. The production of large number of particles 
at z < 1 is a clear indication that we will get here the classical condensate. Let us now discuss the 
role of complex paths in CWKB. 

Case a > 2 I In fig. 5 we plot Nk{z) = N{z)foi five cases. For a = 6.25 the oscillation is not 
transparent because of smallness of Nk{z). We found the oscillations in all the graphs (plotted without 
logscale). We may now easily calculate the period of oscillation setting 6 = w. In terms of z this gives 



-^1/2^2 = expln/yu'^ — almost the same as in [17]. Another aspect of our result is that at small 
z «1 i.e., for long- wavelength mode outside the Hubble radius the occupation number diverges as 



using the expression of Nk{z) for a <2. Thus for physical wavelength Xphy = l/(^?o z) we get 



This result also coincides with Mijic except the factor 1/4. At 2; — > 0, we are to consider the analytic 
continuation as in (39). The emergence of this oscillation also follows from squeezed state formalism 
[4]. We discuss this briefly in concluding section. 

From the discussion of analytical and numerical results, we have the important results that for 
a < 2 we get particle production with large favouring the formation of classical condensate. In 
squeezed state approach to de Sitter spacetime one can show [24], Nk = sinh?rk where is called 
the squeeze parameter. A large is necessary to obtain effective decoherence. 

For a> 2 the modes are well within the horizon as already discussed in this paper just after the eqn 
(52). It has also been shown by Albrccht ct al [4] that in such cases the oscillatory solution becomes 
reasonable approximation. This is what we observe in our CWKB calculation leading to an oscillating 
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FIG. 5. The evolution of one mode occupation number in de Sitter space for different values of mass with 
a > 2. The figure is plotted with cqri (64). The oseillatiou is prominent. For a = 25/4 there are also 
oscillations, not prominent in the figure. Here the particle number is observed to be oscillatory function of 
\ kr] \; there is no formation of classical condensate. 



V. CONCLUSION 

The calculation shows that in de Sitter spacetime the classical condensate can only form for the case 
a < 2 i.e., the inflation field is very weakly coupled and nearly massless. We derive all the essential 
results from a single formula using CWKB modes. We find that the parameter m/Ho plays a vital 
role in the formation of classical condensate . The role of de Sitter horizon distance H^^ and the 
Compton wavelength m = is now clear. Let us generalize the rcsiilt by introducing a nonminimal 
coupling term S,R(fP whose effect is to replace by + 12^iJQ and a by y^ ™ ^^'i^ -^ — 2. For a > 2, 
no classical condensate is formed and for sufficient large z when oscillation dies out the occupation 
number behaves as Planckian spectrum independent of z with 

N = No (e-- - Ij (69) 

where = ^i-f^'^" > i.e., we recover the Planckian spectrum. < is the region for 

the formation of classical condensate. As most of the results obtained in this work substantiate the 
results obtained in [17], without going into further discussion, we let the readers read the concluding 
section of the cited work. What is new is the role of CWKB trajectory in generating the Planckian 
spectrum and discussion of decoherence within the squeezed state formalism. Our CWKB expression 
of particle production is so simple that we need not require to keep track the correctness of the 
expressions employed in each stage for various values of m/Ho [17]. This is the advantage of our 
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approach. Moreover in order to clear up the criticism against the Mijic's works (eqns (26) and (27) 
in the present work), we use Hamiltonian diagonahzation method to obtain the expression of Nk{z) 
vaUd for all z. This needs a minor correction to Mijic's calculation that we have shown without 
changing the basic conclusions of their work. Many feel discomfort with Hamiltonian diagonahzation 
method. Our approach is an example where one can proceed with WKB definition of vacuum instead 
of Hamiltonian diagonahzation. 

A common statement that is made while discussing decolicrcnce is that decoherence is proportional 
to the strength of the coupling and it vanishes if the coupling goes to zero. The analysis carried out 
here for free modes shows that we still get decoherence without being dependent on coupling constant. 
The reason for such a appearance of coupling-independent decoherence may be due to the some specific 
combination of fluctuations and the inflaton zero mode that can assemble to generate a new effective 
field . It is this object which actually rolls down. This interpretation saves the standard picture 
of how metric perturbation are generated during inflation with the spinodial growth of fluctuations 
dominating the late time behaviour of Bardeen variables for superhorizon modes during inflation. This 
view is expressed in [22]. The reason for studying the particle production through CWKB is twofold. 
This allows us to study parametric and spinodial resonances both in FRW and de Sitter spacetime 
under wide range of initial conditions. It should be pointed out that we find parametric resonances in 
chaotic inflationary scenarios for unbroken symmetry (m^ > 0) whereas spinodial instabilities occur 
in new inflation scenario in broken symmetry (m^ < 0). As this is beyond the scope of the present 
work , we refer the reader to see the role of WKB modes in [26] for a study of out of equilibrium fields 
in infiationary dynamics and its relation to density fiuctuations. We like to point out a connecting 
feature of the present work to the standard approach of decoherence in which classicality is studied 
through the squeezed mode parameters. We show that the emergence of squeezed modes is a very 
natural outcome in the CWKB framework with oJk having turning points at real or complex time, 
hitherto not discussed in the literature. 

To relate our work to decoherence, we note that proceeding as in [4], one can easily show that the 
squeeze parameter rfe becomes large because is large at 2; — > for /Hq < 2 since = sin^hrk- 
Not only that we have been able to generate oscillating Nk for rnP' /Hq > 2 (also obtained in [4]) along 
with the emergence of thermal spectrum with correct de Sitter temperature. In most of the earlier 
works dealing with squeeze state formalism (sec [24]) and in [17] emphasis has been put on largeness 
of Tfc to understand decoherence in terms of large particle production. The emergence of oscillating 
Nh and thermal spectrum have not been discussed effectively. In our view the emergence of thermal 
spectrum plus the large particle production are crucial to understand the emergence of classical world. 
In the present work we show that the CWKB approach reproduce these two aspects very clearly. Not 
only that all the conclusions that one gets through squeeze state formalism also follow in this approach 
which we termed as decoherence without decoherence. 
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